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Abstract-  Electrothermal-chemical  (ETC)  ignition 
systems  have  been  demonstrated  in  gun  systems  to 
provide  desirable  characteristics  including  reproducible 
shorter  ignition  delays.  The  optimum  combination  of 
capillary  tube  and  fuse  wire  properties  has  not  been 
identified  yet.  We  present  a  combined  theoretical  and 
experimental  study  of  the  capillary  discharge  with  an  aim 
to  develop  a  capillary  plasma  source  with  efficient 
energy  conversion.  The  major  emphasis  in  the  present 
capillary  discharge  model  is  the  ablation  phenomenon. 
Consideration  is  given  to  different  characteristic  sub- 
regions  near  the  ablated  surface:  namely,  a  space-charge 
sheath,  a  Knudsen  layer,  and  a  hydrodynamic  layer.  A 
kinetic  approach  is  used  to  determine  the  parameters  at 
the  interface  between  the  kinetic  Knudsen  layer  and  the 
hydrodynamic  layer.  Coupling  the  solution  of  the  non¬ 
equilibrium  Knudsen  layer  with  the  hydrodynamic  layer 
provides  a  self-consistent  solution  for  the  ablation  rate. 
According  to  the  model  predictions,  the  peak  electron 
temperature  is  about  1.4  eV,  the  polyethylene  surface 
temperature  is  about  700  K,  and  the  pressure  is  about  10 
MPa  in  the  case  of  a  0.6  kJ  discharge.  In  parallel,  a 
parametric  experimental  study  of  the  capillary  ablation 
process  is  conducted.  The  ablation  rates  are  measured  for 
capillary  tubes  made  of  polyethylene  and  Teflon.  Both 
experimental  measurements  and  simulations  indicate  that 
the  ablated  mass  increases  with  the  peak  discharge 
current  and  that  a  smaller  diameter  capillary  yields  a 
larger  ablated  mass.  It  is  found  that  model  predictions 
agree  well  with  experimental  measurements. 
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I.  INTRODUCTION 


Electrothermal-chemical  (ETC)  guns  are  designed  to  introduce  electrical  energy  into  the 
chamber  of  a  conventional  chemical  gun  system.1  An  ETC  gun  has  the  potential  to 
improve  the  performance  of  conventional  guns  by  enhancing  the  ignition  and 
combustion  characteristics  of  propellants.  One  of  the  significant  advantages  of  the  ETC 
system  is  a  smaller  ignition  delay  time  as  compared  with  conventional  igniters. 
Typically,  ETC  ignition  delay  is  of  the  order  of  1-2  ms  while  advanced  conventional 
igniters  have  been  shown  to  exhibit  6-8  ms  ignition  delay.1  An  ETC  system  has  two 
components,  namely  a  plasma  source  and  a  chemical  propellant  system.  In  this  system, 
the  propellant  is  usually  ignited  with  plasma  introduced  by  a  jet  from  a  capillary 
discharge  or  by  an  arc  from  an  exploding  wire  inside  the  propellant  charge.1’2,3  Several 
plasma  properties  are  important  for  this  application.  The  temperature  of  the  plasma  is 
about  10,000  K.  Significant  radiation  from  the  plasma  leads  to  photochemical  changes 
within  the  propellant  surface  and  in-depth  structure.  The  presence  of  the  various  ion  and 
neutral  species  may  lead  to  chemical  reaction  at  the  propellant  surface.  In  addition,  the 
highly-directed  velocity  of  the  plasma  jet  affects  the  heat  transfer.  All  of  the  above 
plasma-propellant  interactions  are  dependent  on  the  plasma  properties  that  are 
determined  by  the  plasma  source.  Therefore,  accurate  prediction  of  plasma  properties 
from  the  capillary  discharge  is  important  for  understanding  plasma-propellant 
interactions.  The  main  physical  processes  occur  in  the  capillary  in  a  similar  way  to  an 
ablation-controlled  discharge. 

Significant  efforts  have  been  applied  to  model  the  ablation-controlled  discharge  that  is 
commonly  used  as  a  plasma  source  for  ETC.4,5,6  Some  initial  theoretical  models  aimed 
to  analyze  existing  experimental  data  and  build  a  phenomenological  description  of  the 
ablation  processes.  Later,  more  sophisticated  steady  state7,8  and  unsteady  theoretical 
models9,10,11  were  proposed  to  describe  capillary  discharges  for  an  electrothermal  gun. 
However,  these  models  are  based  on  various  assumptions  about  the  ablation  rate.  In 
most  models,  the  ablation  rate  is  determined  somewhat  arbitrarily  or  adjusted  to  fit  the 
measured  total  ablated  mass.  It  was,  however,  shown  recently  that  material  ablation  into 
the  discharge  plasma  is  a  rather  complicated  phenomenon.12,13  For  instance,  the  same 
ablation  rate  can  correspond  to  very  different  discharge  conditions.  A  self-consistent 
model  of  the  capillary  discharge  that  takes  into  account  ablation  kinetics  was 
developed.14,15  The  model  is  based  on  a  kinetic  description  of  the  Knudsen  layer  and  a 
hydrodynamic  description  of  the  collision-dominated  plasma  region.  In  this  paper,  we 
present  a  capillary  model  for  an  ETC  that  includes  self-consistent  consideration  of  the 
ablation  phenomena.  In  addition  we  present  experimental  measurements  of  the  capillary 
ablation  mass  for  various  discharge  conditions  and  capillary  geometry.  We  compare  our 
model  predictions  with  measured  data  in  an  effort  to  validate  the  overall  capillary 
discharge  model. 

II.  EXPERIMENTAL  SET-UP 

Figures  1  and  2  illustrate  the  experimental  fixture  utilized  in  the  present  work.  A  25.4- 
mm  (1-inch)  diameter  outer  insulator  contains  a  6.4-mm  diameter  hole  which  houses  a 
variable  length  capillary  tube  and  inner  insulator.  Both  30-mm  and  38-mm  long 
capillary  tubes  are  investigated.  The  capillary  tube  outer  diameters  are  fixed  at  6.2-mm, 


however,  two  different  inner  diameters  (3.18-mm  and  4.32-mm)  are  considered.  The 
two  different  capillary  tube  lengths  are  accommodated  by  adjusting  the  total  length  of 
the  anode  (brass  rod  +  carbon  disk)  and  the  supporting  inner  insulator  (see  Figs  1,2.). 
The  capillary  nozzle  is  made  of  tantalum/tungsten  and  has  the  following  dimensions: 
length  =  1 1-mm,  OD  =  6.2-mm,  ID  =  4.4-mm.  The  nozzle  serves  as  the  cathode  which 
is  electrically  connected  to  the  anode  by  a  0.10-mm  (0.004-inch)  diameter  aluminum 
fuse  wire.  The  mass  of  the  aluminum  fuse  wire  is  0.6  mg  and  1 .0  mg  for  the  23-mm  and 
38-mm  capillary  tubes,  respectively. 


The  capillary-sustained  plasma  characteristics  observed  are  produced  by  an  electrical 
discharge  directed  across  the  aluminum  fuse  wire.  The  system,  which  delivers  the 
current  across  the  fuse  wire,  includes  a  standard  pulse  forming  network  (PFN)  power 
supply.  The  PFN  power  supply  is  used  to  charge  a  1700  pF  capacitor.  The  system 
charging  voltage  ranges  from  680  V  to  1210  V  with  corresponding  stored  energies  of 
370  J  to  650  J  in  the  experiments  discussed  here.  The  current  is  delivered  across  the 
fuse  wire  producing  the  electrically  charged  plasma.  For  this  system,  the  resulting 
plasma  pulse  width  is  typically  around  385  ps.  The  discharge  current  is  measured  by 
integrating  the  output  of  a  Rigowski  coil  and  high  voltage  probes  are  used  to  measure 
the  voltage  across  the  electrodes.  The  carbon  disk  portion  of  the  anode  is  introduced 
specifically  because  it  has  been  shown  to  have  minimal  erosion  (compared  to  other 
materials)  characteristics  in  the  plasma  environment.  Thus,  the  carbon  disk  has  minimal 
contribution  to  the  nearby  capillary  sustained  plasma.  Once  the  plasma  is  formed  inside 
the  capillary  tube,  the  local  pressure  begins  to  rise  and  the  plasma  vents  through  the 
nozzle  into  open-air.  The  plasma  is  temporarily  sustained  and  enhanced  by  ablating 
material  from  the  capillary  wall.  The  nozzle  mass  is  also  recorded  before  and  after  each 
experiment. 


Figure  2.  Photograph  of  Experimental  Setup 


III.  THE  ABLATION-CONTROLLED  DISCHARGE  MODEL 

In  this  section,  the  model  is  described  for  the  plasma  generation  processes  (ablation, 
heating,  radiation,  ionization  etc.)  and  plasma  acceleration  along  a  capillary  of  a  pulsed 
electrical  discharge.  Figure  3  shows  some  characteristic  regions  in  the  interface  between 
the  discharge  plasma  and  the  dielectric  wall  such  as  an  electrical  sheath  near  the 
dielectric,  the  Knudsen  and  hydrodynamic  layers,  and  a  quasi-neutral  plasma.  Different 
kinetic  and  hydrodynamic  phenomena  determine  the  main  features  of  the  plasma  flow 
including  Joule  heating,  radiative  and  convective  heat  transfer  to  the  dielectric,  and 
electrothermal  acceleration  of  the  plasma  up  to  the  sound  speed  at  the  cavity  exit.  The 


central  region  is  the  quasi-neutral  plasma  that  occupies  almost  the  entire  capillary  since 
typically  the  transition  region  scale  length  is  much  smaller  than  the  capillary  radius.  The 
plasma  region  is  separated  from  the  dielectric  surface  by  the  vapor  layers  (Knudsen 
layer  and  hydrodynamic  layer).  Finally,  the  plasma- wall  transition  region  includes  the 
electrostatic  sheath  attached  to  the  wall. 


Anode  Polyethylene  Cathode  Polyethylene 


Figure  3.  Schematic  of  the  problem  geometry  (not  to  scale)  and  multi-layer  structure  near  the 
ablated  surface. 


The  plasma  is  heated  due  to  electric  current  flowing  through  the  capillary.  The  energy 
transfer  from  the  plasma  column  to  the  capillary  wall  consists  of  the  heat  transfer  by 
particle  fluxes  and  radiation  heat  transfer.  Energy  is  absorbed  by  the  capillary  walls  and 
dissipated  by  thermal  conductivity  and  material  evaporation.  Ablated  flux  builds  up  a 
vapor  layer  in  the  vicinity  of  the  wall.  Under  the  considered  conditions  the  mass, 
momentum  and  energy  conservation  equations  have  the  following  forms: 
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continuum  spectrum  based  on  a  theoretical  model16.  According  to  Ref.  17,  the  radiation 
in  continuum  from  a  carbon  plasma  in  the  considered  parameter  range  provides  the 
main  contribution.  In  addition,  a  black-body  radiation  from  the  plasma  will  be 
considered  and  compared  with  other  energy  sources  and  sinks.  The  particle  convection 
flux  Qf  includes  energy  associated  with  electron  and  ion  fluxes  to  the  dielectric  wall 
that  leads  to  plasma  cooling.  Radiation  and  convection  heat  fluxes  from  the  plasma  to 
the  cavity  wall  (see  Fig.  3)  determine  the  thermal  regime  of  the  capillary  walls.  The 
temperature  inside  the  dielectric  can  be  calculated  from  the  heat  transfer  equation: 

8T  d2T 

—  =  a — y-  (4) 

dt  dr 2 

where  a  is  the  thermal  diffusivity.  Having  calculated  the  plasma  density  and  electron 
temperature  (Eqs.  1-3)  one  can  calculate  the  chemical  plasma  composition  considering 
Local  Thermodynamic  Equilibrium  (LTE)  in  the  way  described  previously18,19.  In  the 
considered  range  of  electron  temperature  (1-2  eV)  and  plasma  density  (1024-1026  m"3)  it 
is  assumed  that  polyatomic  molecules  C2H4  are  fully  dissociated  and  we  will  start  our 


consideration  from  the  point  when  we  have  a  gas  containing  atoms  of  C  and  H.  The 
distribution  of  the  atoms,  multiply  ionized  ions,  and  electrons  obeys  the  corresponding 
Saha  equation.  Knowing  the  plasma  pressure  and  temperature,  ionization  potential,  and 
partition  functions  for  z-fold  ionization  processes  allows  calculation  of  the  atom  and  ion 
species  density.  The  Saha  equations  should  be  supplemented  by  the  conservation  of 
nuclei  and  quasi-neutrality  in  order  to  calculate  plasma  composition.  Plasma  may  depart 
from  the  ideal  state  under  the  considered  high-density  conditions.  Due  to  non-ideal 
plasma  effects,  corrections  are  considered  to  the  ionization  energies,  partition  functions 
and  plasma  pressure.20  The  capillary  wall  ablation  is  modeled  in  the  framework  of  the 
previously  developed  kinetic  model.12,13.  Two  different  layers  between  the  ablated 
surface  and  the  plasma  bulk  are  considered  as  shown  in  Fig.  3:  (i)  a  kinetic  non¬ 
equilibrium  layer  adjacent  to  the  surface  with  a  thickness  of  about  one  mean  free  path; 
and  (ii)  a  collision-dominated  layer  with  thermal  and  ionization  non-equilibrium.  The 
plasma-wall  transition  layer  includes  also  an  electrical  sheath  described  below.  Solving 
the  conservation  equations  in  the  hydrodynamic  layer12,13  the  ablation  rate  has  the 
following  form: 


r=  mVini  =  mnx 


1  (kn2T2  - knxTx) 
m  ,  n ,2 
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It  should  be  noted  that  this  equation  has  a  real  solution  if  the  following  conditions  are 
fulfilled:  k T2fi2  <kTini  and  ni>fi2',  otherwise,  the  solution  of  equation  5  breaks  down. 
Physically,  this  means  that  backflux  is  much  higher  than  the  primary  flux.  Therefore,  in 
this  case,  deposition  on  the  dielectric  surface  takes  place  instead  of  ablation.  The  system 
of  equations  is  closed  if  the  equilibrium  vapor  pressure  can  be  specified  that  determines 
parameters  (ns  and  Ts)  at  the  dielectric  surface.  Various  models  and  an  experiment21 
show  that  the  vapor  pressure  for  polyethylene  can  be  approximated  as  follows: 

P  =  exp(zl-[l-i])  . (6) 

B  1 


where  A=3759;  B=453,  and  T  is  in  K.  Near  the  dielectric  surface  in  the  capillary,  an 
electrostatic  sheath  is  formed  to  maintain  the  floating  potential  on  the  wall  with  respect 
to  the  plasma.  The  potential  drop  of  the  electrostatic  sheath  near  the  dielectric  wall  is 
negative  in  order  to  repel  the  excess  thermal  electrons,  so  that  the  random  electron 
current  density  jeth  is  equal  to  the  ion  current  density  j,. 


IV.  RESULTS 


The  electron  temperature  distribution  during  the  discharge  is  shown  in  Fig.  4.  For 
reference,  the  experimentally  measured  current  waveform  adopted  in  the  simulations  is 
also  shown.  It  can  be  seen  that  the  electron  temperature  temporal  distribution  follows 
the  current  waveform.  The  temperature  peaks  at  1.4  eV  and  then  decreases  towards  the 
end  of  the  pulse. 


Figure  4.  Temporal  variation  of  the  electron  temperature  and  current  waveform  used  in  the 
simulations. 

The  axial  and  temporal  distribution  of  total  pressure  is  shown  in  Fig.  5.  It  can  be  seen 
that  pressure  peaks  at  about  10  MPa  and  then  decreases  towards  the  end  of  the  pulse  to 
about  1  MPa.  Spatial  and  temporal  distributions  of  the  different  species  are  shown  in 
Fig.  6.  All  species  densities  peak  at  about  150-200  ps.  One  can  see  that  the  neutral 
density  is  much  higher  than  the  ion  density,  thus  indicating  that  the  ionization  degree  is 
small  (approximately  10-20  %  for  the  main  portion  of  the  pulse). 
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Figure  5.  Temporal  and  spatial  variation  of  the  pressure  in  the  capillary. 
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Figure  6.  Temporal  and  spatial  variation  of  the  chemical  species  during  discharge  pulse. 


The  surface  temperature  distribution  is  shown  in  Fig.  7.  The  surface  temperature  sharply 
increases  initially  and  peaks  at  about  700  K.  Only  slight  variation  of  the  surface 
temperature  along  the  capillary  is  found  in  accordance  with  previous  predictions. 

Capillary  geometry  determines  the  ablation  process  in  the  capillary  by  affecting  the 
mass  and  energy  balance.  The  dependence  of  the  ablated  mass  on  the  peak  discharge 
current  is  shown  in  Fig.  8.  The  ablated  mass  has  a  linear  dependence  on  the  discharge 
peak  current.  For  comparison,  experimental  data  are  also  shown  in  Fig.  8.  It  can  be  seen 
that  both  simulations  and  experiment  suggest  that  the  ablated  mass  increases  as  the 
capillary  inner  diameter  decreases.  This  effect  is  explained  by  increased  current  density 
in  the  capillary,  which  affects  the  ablation  rate  through  Joule  heating.  Reasonable 
agreement  is  obtained  although  the  predicted  ablated  mass  is  lower  than  that  measured 
in  the  experiments.  It  should  be  pointed  out  that  above  calculations  were  performed 
assuming  that  returned  atoms  and  ions  do  not  form  film  at  the  polyethylene  surface.  On 
the  other  hand  it  can  be  considered  that  only  carbon  atom  and  carbon  ion  deposition 
takes  place,  as  the  hydrogen  atoms  and  ions  will  be  re-evaporated15.  The  last  assumption 
is  supported  by  previous  studies  of  dielectric  (Teflon)  ablation  into  C-F  plasmas 
indicated  that  dielectric  can  be  significantly  carbonized  (charred)  dependent  on 
operational  conditions.22  In  this  paper  we  also  study  parametrically  effect  of  the 
condensation  at  the  dielectric  surface.  We  introduce  a  new  parameter,  v,  which  is  the 


fraction  of  the  backflux  that  condense  at  the  surface.  In  the  case  v=l  all  particles 
returning  to  the  surface  will  condense,  while  small  parameter  v  means  that  only  a 
fraction  of  the  backflux  form  the  film.  Dependence  of  the  ablated  mass  on  this 
parameter  is  shown  in  Fig.9.  One  can  see  the  best  agreement  with  experimental  data  is 
achieved  if  v=0.6-0.7. 
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Figure  7.  Temporal  and  spatial  variation  of  the  polyethylene  surface  temperature. 


Figure  8.  Dependence  of  the  ablated  mass 
(polyethylene)  on  the  peak  discharge  current. 
Comparison  experiment  and  simulations. 
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Figure  9.  Dependence  of  the  ablated  mass 
(polyethylene)  on  the  peak  discharge  current 
with  deposition  fraction  as  a  parameter. 
Comparison  experiment  and  simulations. 


IV.  CONCLUDING  REMARKS 

Self-consistent  analysis  of  the  different  characteristic  regions  near  the  ablated  surface 
allows  determination  of  the  rate  of  ablation  without  employing  various  assumptions 
used  previously.  For  instance,  in  the  present  analysis  it  is  not  necessary  to  assume  the 
energy  fraction  causing  capillary  wall  ablation.  Both  experimental  measurements  and 
simulations  indicated  that  the  ablation  rate  increased  nearly  linearly  with  discharge  peak 


current.  In  addition,  it  was  found  that  the  ablated  mass  increased  with  a  decrease  in 
capillary  radius.  The  effect  of  the  ablated  mass  dependence  on  the  current  density  in  the 
capillary  can  be  explained  by  invoking  the  fact  that  the  main  heat  source  inside  the 
capillary  is  the  Joule  heat  that  depends  on  the  current  density.  Generally,  good 
agreement  was  found  between  model  predictions  and  experimental  measurements. 
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